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> ■ ABSTRACT 

^ : 

^ I We have estimated fluxes of neutrinos and gamma-rays that are generated 

O ! from decays of charged and neutral pions from a pulsar surrounded by supernova 

2 I ejecta in our galaxy, including an effect that has not been taken into consider- 

"^H I ation, that is, interactions between high energy cosmic rays themselves in the 

Qh| nebula flow, assuming that hadronic components are the energetically dominant 

O I species in the pulsar wind. Bulk flow is assumed to be randomized by pass- 

I ing through the termination shock and energy distribution functions of protons 

I and electrons behind the termination shock are assumed to obey the relativistic 

I Maxwellians. We have found that fluxes of neutrinos and gamma-rays depend 

^ ' very sensitively on the wind luminosity, which is assumed to be comparable to 

' the spin-down luminosity. In the case where B = lO^^G and P = 1ms, neutrinos 

should be detected by km^ high-energy neutrino detectors such as AMANDA and 
IceCube. Also, gamma-rays should be detected by Cherenkov telescopes such as 
CANGAROO and H.E.S.S. as well as by gamma-ray satellites such as GLAST. 
On the other hand, in the case where B = lO^^G and P = 5ms, fluxes of neutri- 
nos and gamma-rays will be too low to be detected even by the next-generation 
detectors. However, even in the case where B = lO^^G and P = 5ms, there is 
a possibility that very high fluxes of neutrinos may be realized at early stage of 
a supernova explosion [t < lyr), where the location of the termination shock is 
very near to the pulsar. We also found that there is a possibility that protons 
with energies ~ 10^ GeV in the nebula flow may interact with the photon field 
from surface of the pulsar and produce much pions, which enhances the intensity 
of resulting neutrinos and gamma-rays. 
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Subject headings: acceleration of particles — neutrinos — shock waves — pulsars: 
general — stars: winds, outflows — gamma-rays 

1. Introduction 

It has been about 35 years since Goldreich and Julian (1969) pointed out that a ro- 
tating magnetic neutron star generates huge electric potential differences between different 
parts of its surface and, as a result, should be surrounded with charged plasma, which is 
called a magnetosphere. Gunn and Ostriker (1969) also pointed out the possibility that 
a rotating magnetic neutron star may be a source of high energy cosmic rays. Such high 
energy cosmic rays are considered to be driven along magnetic field lines since these lines 
do not cross. Also, since part of the magnetic fields around the neutron star passes through 
the light cylinder, which means such magnetic fields are open, accelerated charged par- 
ticles are driven to outside of the hght cylinder, which are called as pulsar winds. The 
pulsar winds are usually considered to be composed of electron-positron pairs since electron- 
positron pairs will be created so as to eliminate electric fields that are parallel to mag- 
netic fields at a region where net charge density is not equal to the Goldreich-Julian den- 
sity (Ruderman and Sutherland 1975; Shibata 1991). However, it is also pointed out that 
hadronic component may exist in pulsar winds as a consequence of the net charge neutrality 
in the outfiow (Hoshino et al. 1992; Bednarek and Protheroe 1997; Protheroe et al. 1998; 
Blasi et al. 2000; Amato et al. 2003). Moreover, it is pointed out that hadronic components 
may be the energetically dominant species although they are dominated by electron-positron 
pairs in number (Hoshino et al. 1992). This is because inertial masses of hadrons are much 
larger than that of electron. 

Based on the assumption that hadronic components are not negligible in pulsar winds, 
some scenarios are proposed to produce high energy neutrinos and gamma-rays from decays 
of charged and neutral pions that are produced by interactions between hadronic, acceler- 
ated high energy cosmic rays and surrounding photon fields and/or matter. Atoyan and 
Aharonian (1996) estimated fiux of gamma-rays and discussed its contribution on the ob- 
served spectnim of the Grab nebula, although they concluded that its contribution may be 
important only at energies above 10 TeV. Bednarek and Protheroe (1997) proposed that ac- 
celerated heavy nuclei can be photo-disintegrated in the pulsar's outer gap, injecting energetic 
neutrons which decay into protons. The protons from neutron decay inside the supernova 
ejecta should accumulate, producing neutrinos and gamma-rays in collisions with the matter 
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in the supernova ejecta.^ There are some papers based on this scenario and flux of neutri- 
nos and/or gamma-rays is estimated (Protheroe et al. 1998; Bednarek and Bartosik 2003; 
Bednarek 2003a; Bednarek 2003b; Amato et al. 2003). Beal and Bednarek (2002) proposed 
that accelerated cosmic rays will interact with the photon fields inside the supernova remnant 
at the very early phase (within ~1 yr) of a supernova explosion. Bednarek (2001) calculated 
the extragalactic neutrino background based on this scenario. 

In this study, we estimate fluxes of neutrinos and gamma-rays including an effect that 
has not been taken into consideration, that is, interactions between high energy cosmic rays 
themselves. This picture is based on the works given by Rces and Gunn (1974) and Kennel 
and Coroniti (1984). Rees and Gunn (1974) pointed out that the supersonic pulsar wind 
would terminate in a standing reverse shock located at a distance from the pulsar. Beyond 
the termination shock, the highly relativistic, supersonic flow is randomized and bulk speed 
becomes subsonic, obeying the Rankine-Hugoniot relations. According to Hoshino et al. 
(1992), who studied the theoretical properties of relativistic, transverse, magnetosonic col- 
lisionless shock waves in electron-positron-hcavy ion plasmas, proton distribution functions 
in the down stream are found to be almost exactly described by relativistic Maxwellians 
with temperatures Tp^2/li^p 

0.34, where Tp^2 is temperature at the down stream, 71 is 
the bulk Lorenz factor of protons in the up stream. It is noted that protons are not ther- 
malized through the interactions with protons themselves, and/or electrons but just obeys 
the Maxwellian distribution through transferring cyclotron waves. This subsonic flow speed 
would be, by communicating with the nebula boundary at r^v via sound waves, adjusted 
to match the expansion speed of the supernova remnant (that is, supernova ejecta) at the 
innermost region. This subsonic flow is called as nebula flow in the study of Kennel and 
Coroniti (1984). We also adopt this definition in this study. In this study, we calculate flux 
of neutrinos and gamma-rays from charged and neutral pion decays in the nebula flow which 
arc produced through the interactions between high energy protons themselves, assuming 
that energy distribution functions of protons obey the relativistic Maxwellians. 

In this study, as the previous works, we assume that protons are energetically dominant 
in the pulsar winds. Thus, we describe the nebula flow using the proton mass as a unit of 
mass. We estimate pion production rates due to the interactions between high energy protons 
using proper Lorenz transformations. As for the cross sections between protons, scaling 
model (Badhwar et al. 1977) is adopted. Isobar model is not included in this study since 
we consider production rates of high energy pions. Calculating the spectrum of neutrinos 



They call this region as a nebula, which is slightly different from the definition of a nebula given by 
Kennel and Coroniti (1984). In this study, as explained below, we adopt the definition presented by Kennel 
and Coroniti (1984). 
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and gamma-rays in the termination shock rest frame, we estimate these fluxes at the earth 
assuming that the pulsar is located at 10 kpc away from the earth. 

In section 2, we explain the formulation in this study. Results are shown in section 3. 
Discussions are given in section 4. Summary and conclusion are presented in section 5. 



2. FORMULATION 

2.1. Confinement of Pulsar Wind by the Surrounding Supernova Remnant 

In this subsection, we briefly review the steady state, spherically symmetric, magneto- 
hydrodynamic model of nebula flow presented by Kennel and Coroniti (1984), in which a 
highly relativistic pulsar wind is terminated by a strong MHD shock that decelerates the 
flow and increases its pressure to match the boundary conditions imposed by the surround- 
ing supernova remnant. In this study, we consider the possibility that the protons are the 
energetically dominant species in the flow as pointed by Hoshino et al. (1992). According 
to the work presented by Hoshino et al. (1992), the proton distribution functions behind 
the termination shock are almost exactly described by relativistic Maxwellians. Thus we 
can deflne temperature behind the shock wave and solve the following nebula flow using the 
magnetohydrodynamic equations. 



2.1.1. Super-relativistic Magnetohydrodynamic Shock 

It is usually introduced the parameter a, which is the ratio of the magnetic plus electric 
flux to the particle energy flux at just ahead of the termination shock, 

- 

AnnuVmpC^ ' 

where B is the amplitude of the magnetic field in the observer's frame, n is the proper density, 
u is the radial four speed of the flow, F is the bulk Lorenz factor of the flow (F^ = 1 + m^), 
TRp is the proton mass, c is the speed of light. It is assumed that the energy density of 
electric flled is nearly equal to that of magnetic fleld. The luminosity L at just ahead of the 
termination shock is described as 

L = A'KnTurlmpC^{l + a), (2) 

where is the radius of the termination shock. The termination shock is assumed to be 
stationary relative to the pulsar and the observer. Thus the Rankine-Hugoniot relations 
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for 90° shocks can be used directly to obtain the physical quantum in the observer's frame 
behind the shock wave: 



niUi = n2U2 



(3) 



UiBi _ U2B2 



(4) 



EBi EB2 
TiAfi + = r2//2 + 



AnriiUi 



47rn2M2 



(5) 



HiUi H h 



Bl 



P2 

IJ.2U2 H h 



Bi 



(6) 



riiUi SnriiUi n2U2 87^2^2 

Subscripts 1 and 2 label upstream and downstream parameters. E denotes the shock frame 
electric fields and ^ — {e -\- P)/n is the specific enthalpy, which for a gas with an adiabatic 
index 7 is defined by 

n 



jji = mc + 



7 



7-1 



(7) 



In this study, we set 72 to be 4/3, as explained in Appendix A. 

Kennel and Coroniti (1984) introduced a parameter Y to obtain the downstream radial 
four speed, M2. Definition of Y is 

B2 T2U1 



Y = 

Bl V1U2 

Using this parameter and 72 =4/3, Eq. (6) can be expressed as 



(8) 



y2 _Y 



T2U2 
27rnimpC^ Ui 



2 , A ^1 
' 4; Ti 



+ 



r2«2 v''' + 4 



Bl r\ 



87r(ni/XiMf + Pi 



Although there are some trivial typos for the expression of Eq. (9) in the original paper, 
they does not affect the following discussion. 

Assuming that the bulk Lorenz factor of the pulsar wind is sufficiently large and the flow 
is sufficiently cold, we can approximate Fi, Pi, and /xi as Fi ~ ui, Pi ~ 0, and //i ~ m^c^. 
Using these approximations, we can rewrite Eq. (9) as 

2 / 1 \ 2 



2 / 2 
^2 ^2 + 



[I + ul) ul - 



1 



41 + (7 



(10) 
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Solving Eq. (10), we can obtain the radial four speed behind the shock wave as 

o 8a^ + 10(7+1 1 



1/2 



(11) 



, , + — ^ \64:a\a + 1)' + 20a(a + 1) + l] 

16(c7+ 1) 16((7+ 1) L V y V y J 

Note that the radial four velocity behind the shock wave is determined by only one parameter, 



a. 



As for the pressure P2, it can be obtained as 

P2 1 



where we assumed that 112 ^ 4(P2/'^2)- 



l + a|l-^^ 

U2 



(12) 



In this study, we adopt a slightly different way to determine the temperature from the 
way presented by Kennel and Coroniti (1984), using the results of Hoshino et al. (1992). Ac- 
cording to Hoshino ct al. (1992), the distribution functions of protons behind the termination 
shock are almost exactly described by relativistic Maxwellians as 



N{'y) = A'jexp 



■^(7-1) 



(13) 



where 7 is the Lorenz factor that represents the random motion (note that this is not the 
adiabatic index), /c^ is the Boltzmann constant and A is the normalization factor in units of 
cm~^. In this case, the relation 



P2 = ^n2kBT2 



(14) 



holds exactly (see Appendix A in detail). Thus the temperature T2 can be derived from 
Eqs. (12) and (14) as 



sr. 



l + a(l-^^ 

U2 



(15) 



As a result, we can estimate the energy density (e) [erg cm ^] as 

/oo 
72e-"^d7, (16) 

where a — rripC? j ksT. Since we can obtain the physical quantum just behind the termination 
shock, we consider the resulting flow from the termination shock to the remnant. This flow 
is usually called as nebula flow. 
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2.1.2. Nebula Flow 

The steady and spherically symmetric nebula flow is described by the equations men- 
tioned below. These are; the equation describing the conservation of number flux, 

^icnur^) = 0; (17) 
dr 

one describing the conservation of magnetic flux in the magnetohydrodynamic approxima- 
tion, 

^ 0; (18) 



dr V r 

one describing the propagation of thermal energy. 



^iur^e)+P^{r''u)^Q- (19) 
dr ^ ' dr 



and one describing the conservation of total energy. 



d_ 
dr 



0. (20) 



Note that the notation for the thermal energy density is shghtly different from the original 
paper. 

When the relation e = 3P holds, Eq.(19) reduces to 

:^lnf-^Uo. (21) 



dr \ n'^/^ ^ 

These are the basic equations describing the nebula flow which connects the condition behind 
the termination shock and the inner-edge of the supernova remnant. 

Here we assumed that the distribution functions of protons remain Maxwellian and 
temperature can be determined at each radius of the nebula flow. However, we should adopt 
the Boltzmann equations coupled with the evolution of magnetic fields to obtain the exact 
solution for the nebula fiow, since it is not apparent for protons to transfer their energies 
with each other enough to obey the Maxwellian distribution. However, as shown in section 3, 
charged and neutral pions are mainly produced just behind the termination shock (see also 
figures 8 and 12), where the energy distribution of protons obey relativistic MaxweUian. The 
nebula fiow can be described approximately as a free-expansion and temperature and density 
in the nebula fiow becomes smaller along with radius, which means pions are less produced 
along with radius. Thus we consider that the order-estimate of the neutrino fiux produced 
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by pp collisions will be valid even if the MHD equations are adopted. We will estimate the 
neutrino flux using the Boltzmann equations coupled with the evolution of magnetic fields 
in the forthcoming paper. 

The solution for the nebula fiow can be obtained analytically (see Appendix C for 
details). Using the solution, total pressure Pt in the postshock can be expressed as 



vzY'^' + 4(1 + 



Sn 47rr^c(l + a) 
where v is defined as v = u/u2 and z is defined as 2; = r/r 



2uW 



(22) 



As explained in section 2.1.3, wc investigate the case a <^ 1, because the speed of the 
nebula flow at the boundary between the remnant and nebula flow is set to be almost same 
with the speed of the remnant (~ 2000 km s~^). In the case of cr <C 1, 



U2 



l + 9a 
8 



1/2 



and 



(23) 



';{l-7a). 



Thus the total pressure can be written approximately 



■{niTiUimpC^) 



+ ^ 



1 + 



(24) 



(25) 



This approximation is used in section 2.1.3. 



2.1.3. Boundary Conditions 

As for the inner boundary condition, the parameters are bulk Lorenz factor and lu- 
minosity of the wind. As for the bulk Lorenz factor, we can estimate its upper limit as a 
function of the rotation period and the amplitude of the polar magnetic field of the pul- 
sar (Goldreich and Julian 1969). The corotating magnetosphere is bounded by a field fine 
whose feet are at sin^o ~ (fi-R/c)^/^ = {2R/Py/^ x 10~^, where 9o is the zenith angle (in a 
unit of rad), Q is the angular velocity of the pulsar, R is the radius (in a unit of cm), and P 
is the rotation period (in a unit of s) . The potential difference between 9o and the pole is 

A* = i(^)'i?i^p (26) 
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for c <^ 1 when the surrounding region of the pulsar is vacuum. Here Bp is the amphtude 
of the polar magnetic field in units of G. Although the magnetosphere should be filled with 
plasma whose number density is the Goldreich- Julian value, the most energetic escaping 
particles can be estimated using Eq. (26). It will be ~ A$/2, or 



where Z is the atomic charge. Thus the bulk Lorcnz factor of protons can be estimated as 
r„„^3.2xlO»fJ^Ul^V7^V. (28) 



W^G j \ P j VlO^cm 

In this study, we assume that the bulk Lorenz factor of the pulsar wind is monolithic and 
its upper limit is given by Eq. (28). 

The wind luminosity can be estimated when it is assumed to be comparable to the 
pulsar's spin down luminosity. Under the assumption, it can be expressed as 

^43. ( \ ( \ ( ^''^^ \ -1 



i = 9.6xl0'^^^j ^I^j ergs-. (29) 

In this case, the angular velocity evolves as a function of time as 

/ R2 P6q2 \ -1/2 

Qit) ^nAl + -^t rad s-\ (30) 

\ 3c''/ / 

where f2j and / ~ 10"^^ g cm^ arc initial angular velocity and inertial moment of a pulsar, 
respectively. We can estimate the spindown age which is defined as the time when the 
angular velocity Vt{t) becomes 20.^ as 

Win^6.2xl0M^ [——] — - s. (31) 



Bp ) \ R J Vims 

In figure 1, we show the spin down age [yr] as a function of period of a pulsar (solid line) as 
well as the spin down luminosity (dashed line) . Amplitude of the magnetic field at the pole 
of a pulsar is set to be lO^^G in this study. 

Next, we consider the outer boundary condition. It is naturally considered that the 
inner region of a supernova remnant is composed of He layer that includes heavy nu- 
clei (Hashimoto 1995; Woosley and Weaver 1995; Thielemann et al. 1996; Nagataki et al. 1997; 
Nagataki 2000) whose escape velocity is ~ 2000 km s~^ for the free-expansion phase (Haas et al. 1990; 
Spyromilio et al. 1990). In this phase, the initial explosion energy is almost all manifested 
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by kinetic energy of the ejected matter; thermal energy comes to 2 or 3% of the initial ex- 
plosion energy. Thus we model the innermost region of the ejecta as follows. Wc consider 
the 6Mq He layer which is the typical mass of helium for the progenitor of collapse-driven 
supernova (Hashimoto 1995). Then, the speed of the escaping velocities of the He layer is 
2000 km s~^ (=Knin) for the inner edge and 3000 km s~^ (=Knax) for the outer edge. The 
thermal energy in the He layer is assumed to be 

E,^ = 0.02 X 10^1 X erg 
* 20Mq ^ 

= 6 X 10^^ erg, (32) 

where 0.02 represents the fraction of the thermal energy in the ejecta, 10^^ is the typical 
explosion energy of a collapse-driven supernova, 6Mq is the mass of the He layer and 2OM0 
is the typical total mass of the ejecta. Thus the volume occupied by the He layer can be 
calculated as 



' 2" I' max minj \ 



t 



Isec 



(33) 



where t is the age of the supernova remnant. Thus the typical pressure in the He layer can 
be estimated by solving the equations 

^th = ^(^e + NHe)kBT + SaT^V (34) 

and 

P = (rie + nue)kBT + aT\ (35) 

where A'^c and A^hc ai'e the total number of electron and helium in the He layer, ric and 
are the number density of electron and helium, and a is the radiation constant, respectively. 
Here we approximated that the He layer is composed of helium, neglecting the contamination 
of heavy nuclei. Using this model, we can estimate the pressure at the innermost region as a 
function of time and the result is shown in figure 2. The discontinuity at i ~ 5 yr reflects the 
transition from photon-dominated phase to matter-dominated phase. This happens when 
the optical depth of the supernova ejecta becomes lower than unity and pressure of photon 
fields is set to be zero. 

As for the Crab nebula (age of the Crab is ~ lOOOyr), the pressure at the innermost re- 
gion of the remnant is estimated to be in the range (1-10) x 10^ dyn cm~^ (Kennel and Coroniti 1984) . 
On the other hand, Iwamoto et al. (1997) estimated the pressure of the innermost region of 
the ejecta numerically and reported that P ~ 10^^ dyn cm~^ a,t t — 100 sec (~ 3 x 10~^yr). 
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We can find from figure 2 that the estimation by Eq. (35) reproduces these values fairly well. 
Thus we adopt this formula in this paper throughout. 

Finally, we consider how to determine the position of the termination shock. The termi- 
nation shock should be initiated at the innermost region of the remnant where the relativistic 
pulsar wind hits. Then, if the stationary state exists, the position of the termination shock 
propagates inward so as to attain the pressure balance between the nebula flow and rem- 
nant. If the pressure of the nebular flow (P„) can not be comparable to that of the remnant 
(Pr), what happens? If (a) P„ > Pr, the nebula flow should push the remnant until the 
pressure balance is attained. In such a situation, interactions between relativistic protons 
in nebula flow and non-relativistic protons in the ambient remnant should be effective and 
may result in effective production of high energy neutrinos. This situation will be similar to 
that of the previous works such as Bednarek and Protheroe (1997). In such a situation, the 
interactions between the relativistic protons in the nebula flow and photons in the remnant 
may also become important, which will be similar to the situation that Beall and Bednarek 
(2002) presented. Such a situation is outscope of this study and detailed estimation of flux 
of high-energy neutrinos has been done by the previous works. On the other hand, if (b) 
P„ < Pr, the reverse shock initiated at the inner edge of the remnant would be driven back 
to the pulsar, as pointed out by Kennel and Coroniti (1984). This situation is within our 
scope and considered. These situations can be distinguished by using Eq. (25). The total 
pressure normalized by 2{niriUimpC^)/3 ~ P2 can be solved as a monotonic function of z, 
as shown in figure 3. In this figure, a is set to be 6.7 x 10~^ = (Tc, which is used throughout 
this paper as mentioned below. Prom the definition, z^m — 1 and ^^max = (-Rrem/'^s), where 
Prem = 2 X lO^t IS the radius of the inner edge of the remnant at time = t sec. If the pressure 
in the remnant (Pr) is in the range Pt(-2 = ^max) < Pr ^ Pt{z = z„^m), the radius of the 
termination shock is determined at {R < Vg < Prem)- Otherwise, the situation is (a) 
Pn > Pr or (b) Pn < Pr. 

As for the velocity, the radial four velocity behind the shock wave depends on only one 
parameter, a (Eq. (11)). Also, in the small a limit, the asymptotic solution for 
be written as 

.2 \ 1/2 



^00= (^f^) , (36) 



which corresponds to 



/Soo = T^. (37) 

Thus (7 should be ~ 6.67 x 10~^ in order to obtain the stationary nebula fiow in which the 
velocity of the nebula fiow is almost same with that of the innermost region of the remnant. 
In this study, we set a — ac throughout in this study for simplicity. 



Of course, it is noted that the value of a should be not determined by the outer boundary 
condition, but by pulsar's condition in principle. Thus our requirement, which is also adopted 
in Kennel and Coroniti (1984), may not be so realistic. Thus, before we go further, we 
should discuss here what will happen if (i) o" > o"c or (ii) a < (7^. In the case of (i), the 
speed of the nebula flow is faster than that of the remnant and additional pressure (ram 
pressure) is given to the nebula flow. Thus discussion on the pressure balance mentioned 
above should be modified by introducing the effect of the ram pressure when the ram pressure 
is comparable with or higher than the total pressure. In the case of (ii) , as long as the contact 
discontinuity is adopted as a boundary condition, the location of the shock can not be far 
away from the remnant so that the nebula flow can catch up with the remnant. In this case, 
the resulting neutrino flux in the nebula flow should be small because the volume of the 
cavity (i.e. region between the shock and remnant) should be kept small. Also, there is a 
possibility that shocked region and remnant continue to separate from each other. If so, the 
rarefaction wave will be generated from the remnant and the rarefaction wave will connect 
the remnant and shocked region, achieving the boundary condition of contact discontinuity 
between rarefaction wave and shocked region. It will be important to investigate what kind 
of boundary condition is the most proper as a next step of this study. 



2.2. Microphysics of Proton-Proton Interaction 

In this section, the formulation of production of neutrinos and gamma-rays from pp 
collisions is presented. In section 2.2.1, we calculate at flrst the pion spectrum in the rest 
frame of the nebula flow. Then, the spectrum is transfered to the observer's frame. In 
section 2.2.2, microphysics of pp interactions is explained. Finally, we calculate the flux of 
neutrinos and gamma-rays as products of pion decays in section 2.2.3. 



2.2.1. Formulation of Pion Production 

First, we consider the number spectrum [particles cm~^ s~^ erg~^] of charged and neutral 
pions in the rest frame of the nebula flow. It is calculated by considering two protons with four 
momenta pi and p2, moving towards each other with relative velocity Vj-ei- The number of col- 
lisions that occur in a volume dV, for a time dt, is a frame invariant quantity, which in an arbi- 
trary reference frame can be written as (Landau and Lifshitz 1975; Mahadevan et al. 1997) 



dRi2 = appV.,ei——nin2dVdt 
hii rj2 
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= ccTppmns V^(A - 02)^ - {01 X 02?dVdt, (38) 

where app is the total cross section, rii and n2 are the number density of protons with their 
four momenta are pi and p2-, respectively. Thus the number spectrum of pions [particles 
cm~^ s~^ erg~^] is calculated as 

nE.) _ f^,^^ r%,, r......^'^Pp(7i,72,cos^) 



dV 



r . r , f J ^ c?M7i,72,cosg 

= cy dr/2 I dr/i j dcosO — n(i?, 7i)n(i?, 72) 



X V(/^l-/^2)2-(/^lX/^2)^ (39) 



R 



where 71, 72 are the respective Lorenz factors of the two protons, cos — Pi ■ P2/ 
is the radius with respect to the neutron star, n{R, 7) is the differential number density of 
protons at position R, and dapp{^i, ^2, cos 9) / dE^ is the differential cross section of proton- 
proton interaction, which is explained in section 2.2.2. It is noted that the more energetic 
proton is labeled as 2 and the less energetic one is labeled as 1 in this formulation, which 
means 72 is always larger than 71. In this frame, the distribution of pions in momentum 
space is isotropic. Thus the number spectrum of pions in unit solid angle [particles s^^ crg^^ 
sr~^] is simply expressed as 

dn 47rAy dV ^ 

where Ay is the fluid element. 

Next, we transfer the obtained number spectrum in the fluid- rest frame to the one in 
the observer's frame. It is noted that we consider a receiver in the observer's frame. Thus we 
consider both special relativistic effect and Dopplcr effect. The relative four velocity between 
the nebula flow and observer is considered to be u{r) (i.e. we neglect the effect of proper 
motion of the pulsar with respect to the earth). The result is expressed as (see Appendix B 
for derivation) 

dn r2(l -/3cos^)2 47r ' ^ ' 

where bars are labeled for the quantum in observer's frame, 9 is the angle between the three 
dimensional velocity of the fluid element and direction of the observer measured from the 
fluid element, and f is the bulk Lorenz factor of the fluid element in the observer's frame. 
The relation between E^^ and E^^ is expressed as 

En = -7 i =rK- (42) 

r(l-/3cose) ^ ^ 
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This formulation is also used to calculate the flux of neutrinos and gamma-rays in the 
observer's frame. In this study, however, this transformation is not so effective except for 
the region just behind the termination shock, because the speed of the flow behind the shock 
is not so high in the case of cr <C 1 (see Eqs. (23) and (36)). 



We show the method of calculating the differential cross section of proton-proton in- 
teraction, dapp{ji. j2, cos 9)/ dEjr J which is introduced in section 2.2.1. Wc estimate this 
quantum by two steps. First, we transform the fluid-rest frame to the rest frame of parti- 
cle 1 (i.e. laboratory frame), because formulation of the differential cross section for pion 
production in this frame is presented by a number of previous works (Badhwar et al. 1977; 
Stephens and Badhwar 1981; Dermer 1986a). Using this formulation, we can estimate the 
probability of producing pions with four momentum p'^ in this frame. Next, we transform 
again the coordinate to the fluid-rest frame and obtain the differential cross section of pion 
production in the ffuid-rest frame. 

Let us begin with transforming the ffuid-rest frame to the rest frame of particle 1. Wc 
choose X — y plane so that the particles 1 and 2 moves in this plane. We also choose x axis 
to be aligned with the direction of velocity of particle 1. In fluid-rest frame, the particle 2 
collide with particle 1 with angle 9. In the rest frame of particle 1, this angle becomes 9' . 
Pions are produced as a result of this collision with zenith angle a and azimuthal angle (f)' . 
Here we choose a = axis to be aligned with the direction of velocity of particle 2 in the rest 
frame of the particle 1. We also set = and tt when the pion moves in the x' — y plane. 
We show a sketch of the geometry concerning the individual scattering events in Figure 4. 
The left panel shows the geometry in the fluid-rest frame, while the right panel shows the 
geometry in the rest frame of particle 1. 

In the fluid-rest frame, the four momentum of particle 2 can be expressed as 



2.2.2. Differential Cross Section of Proton- Proton Interaction 
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This is transformed to 
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(44) 



Thus, the relation between d and d can be expressed as 

^2y _ \f' 



tan^ = -- ^ = 

P2x 



'7i 



Isin^ 



71 (-/5i72 + a/72 - 1 cos 6* j 



(45) 



In this frame, the differential cross section for pion production can be expressed as 
(Naito and Takahara 1994) 



dE' 



^jE'^-mld' 



27r 



(icoso; 



(46) 



^dp'^c\ 

where cosa|^^^ = [7^^^ - ^max(s)] / [f^dcP^c] (-1 < cosa^^^ < 1), 7c is the Lorenz fac- 
tor of the center-of-mass system with respect to the rest frame of particle 1, E^^{s) = 
{s — ArripC'^ + m^c^)/2s^/^ is the energy of pion in the center-of-mass system when the pion 
obtains the maximum angle a^ax ^^^^ frame of particle 1 (Dermer 1986b), and s is 

the Mandelstam variables for s-channel. The estimation of the Lorenz invariant quantum 
E'^d^Gpp/ dp'^c^ is explained below. 



The four momentum of produced pion in the rest frame of particle 1 can be expressed 



as 



( 



\ 



\ 



p^{cos a cos 9' — sin a cos (f) sin 9 } 
— p^{cos a sin 9' -\- sin d cos 0' cos 9 } 
sin a sin (f) 



(47) 



It can be transformed to the one in the fluid-rest frame as 

/ 7i 7i/3i 0\ / 7i ^Jp'^ + mlc^ + 7i/5ij9^{cos a cos 9' — sin d cos (f) sin 9'}^ 

7i/3i \/pl^ ~V m\c^ + 7ip^{cos a' cos ^' — sin d cos 0' sin ^'} 
— p^{cos a' sin ^' -|- sin a. cos 0' cos 
Y U U U 1 y y sin a' sin 0' y 
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Thus, since 



J dE^ 



dE^ 



' dapp{E^, E2) 



dEi 



(49) 



dcTppiji, ^2, cos 9) /dE^^ can be expressed as 



(^c^pp(7i)72,cos^) _ d 



dE^ 



dE^ 




dE'^d(f)'d cos a \J E'^ - m^d- ( -E^-rr 



,/ do. 



VP 



' dp'^c^ 



xif(7i-\/p^2 + m^c^ + 7i/3ip^{cos a cos O' — sin a cos 0' sin 9 } — E^) 
xH{Et^ + dE^r — ^w/v^ + m^c^ + 7i/3ip^{cos a' cos ^' — sin a cos 0' sin ^'})}, 



(50) 



where iiffx) is the Heaviside function. 



Here we explain the formulation to estimate the Lorenz invariant cross section E'^^a^j d-p^e" . 
In the center-of-mass system (labeled as *). the invariant cross section is inferred to have a 
form (Badhwar et al. 1977; Naito and Takahara 1994) 



^d^p* (l + 4m2cVs)^ 



1^1 — xY exp 



1 + 4m2cVs 



(51) 



where x 



{xf + (4/s)(pf + m2c'^)}V2 and g = (Ci + C2PI + Cj,f^. Definition of x| is 
Here we decompose the momentum of pion in the center-of-mass 



rripC'^. 



xjj = P|jc/y^s/4 - 

system of protons into the component which is parallel to the velocity of the protons {p*) 
and the one which is perpendicular to the velocity {p*j_). The parameters A, 5, Ci, C2, C3 
and r for neutral and charged pions are tabulated in table 1. Note that the Lorenz invariant 
cross section is expressed in units of mb/ (GeV^/c^), momentum of pion is expressed in units 
of GeV/c, and energy is expressed in units of GeV. This fitting formula is derived from 
the study of pp collisions for incident proton energies 6-1500 GeV. Thus we have to note 
that this formula has to be extrapolated to the high energy range in this study. However, 
as shown in Naito and Takahara (1994), even if the scahng model is extrapolate to higher 
energies, they confirmed that the scaling assumption does not much affect the resulting 
gamma-ray spectrum by comparing their results on the resulting gamma-ray spectrum up 
to ~ 10^ GeV with the gamma-ray production by the mini-jet model including the QCD 
effect (Gaisser and Halzen 1987; Bcrczinsky et al. 1993). As for the low energy range (< 5 
GeV), the isobaric model (Stecker 1971), in which a single pion is produced in a collision, is 
considered to be better than the scaling model. However, we consider much more energetic 
protons and the difference between the isobaric model and the scaling model should be 
unimportant in this study. 
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2.2.3. Flux of Neutrino and Gamma-ray 



The neutral and charged pions decay into gamma-rays, electrons (positrons), and neu- 
trinos as follows: 



TT ' 

7r~ 



TT ^ 7l 72 

» i/^ ^ 6+ i/g + + 

+ J^H ^ e~ + Pe + J^ix + J^n- 



(52) 
(53) 
(54) 



In the case of 2-body decay (tt" — > 71 -|- 72 and tt^ — >• /^^ -|- t/^(i/^)), the four momentum 
of gamma-ray and neutrino in the rest frame of the pion can be obtained easily by calculating 
the conservation of four momentum. As for the 3-body decay (//^ — > -|- t'e(^'e) + ^ni^n)), 
conservation of four momentum and angular momentum (i.e. spin) have to be taken into 
account to estimate the resulting energy spectrum of secondary neutrinos in the rest frame 
of the charged pion. We adopt the formulation presented by Dermer (1986a). It is noted 
that the distribution of neutrino in the momentum space is, after all, isotropic in the pion 
rest frame. 

Next, we transform the obtained number spectrum of neutrino in the rest frame of pion 
to the one in the fluid-rest frame. We can transform the number spectrum of neutrinos in 
pion-rest frame to the one in fluid-rest frame by two steps, that is, Lorenz boost and rotation. 
Let {0,(j)) be the zenith angle and azimuthal angle of pion's velocity in the fluid-rest frame. 
We set xyz coordinate in the fluid-rest frame and x"y" z" in the pion-rest frame. Here we 
choose x" axis to be parallel to the direction of pion's velocity. When the four momentum 
of neutrino in the pion-rest frame is {pi, pi cos a , pi sin a" cos p", pi sin a" sin p"), it is 
transformed as 
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Here a" and /3" are the zenith angle and azimuthal angle of neutrino in the pion-rest frame. 
It is noted that the zenith angle a is defined as the angle between the velocity of ncTitrino 
and x " angle. 7 is the Lorenz factor of pion in the fluid- rest frame. From Eq. (55), the 
energy of neutrino in the fluid-rest frame is expressed as 

-Bj./c = 7p|^(l -F/3C0SQ;"). (56) 

Thus, by considering the conservation of number of neutrino in the fluid element, the relation 

dN, 



(55) 



'Jk^^- 1 1 1 1 1 ^^'^ 



„ dn" (pi) dp" sin a" da , ^ 2^-^^'' ,^ rr,-, ^ "\ 

^ "^^ -dpAT^pl-—dEJ{E^ - 7p^(l + pcosa )c}(57) 
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can be derived. Here dN^/dEi, is the number spectrum of neutrino [particles GeV"^] in 
the fluid-rest frame, dn' {pi) / dpi is the momentum distribution of a neutrino in the pion- 
rest frame [(GcV/c)~^], and d^N'^/dp^ is the number spectrum of pions in the momentum 
space in the fluid-rest frame [particles (GeV/c)~^]. Finally, the number spectrum of neutrino 
[particles GeV~^] in the fluid element in the fluid-rest frame can be expressed as 

' f f f dpl^^ sin a" da" dp/-^6 |l - ^^-^^ +f cosc.")c 1 ^^^^ 



dE, 2E,J J J dp'l dp^ { E, 



It is noted that the number spectrum of pions in the momentum space in the fluid-rest frame 
is isotropic so that it can be reduced to the number spectrum of pion dN'^ /dpT^ [particles 
(GeV/c)-^]. 

Finally, we transform it to the one in the observer's frame, using the formulation pre- 
sented in section 2.2.1. 



2.3. Detection of Neutrinos 

We briefly review the detectability of high energy neutrinos at km^ detector such as 
IceCube. 

In this study, we adopt the formulation presented by Gaisser and Grillo (1987). The 
event rate [events s^^] of the neutrino-induced signal whose energy is larger than E^ at the 
detector with effective area A^s is estimated as 

dN 

S{> E^) = Aeff / dE,--^P{E,, E^), (59) 

where dNi^/dE^, is the neutrino energy spectrum [particles cm^^ s"^ erg~^] and P{E^, E^) is 
the probability that a neutrino aimed at the detector gives a muon with energy above E^ at 
the detector. The latter is given by 

pE^ pE^ 7 ^oo 

P{E^,E,)^ dE; dE^-—NA dXgiX,E;,E^), (60) 
Jb^ Je* Jo 

where A^^ is the Avogadro's number and g(X, E*, E**) is the probability that a muon pro- 
duced with energy E** travels a distance X g cm^^ and ends up with energy in [E*, E*+dE*]. 
Gaisser and Grillo (1987) assumed that g{X, E*, E**) has a form of 



g{X,E:,E:*)^ Xo = -ln 



(61) 
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where a = 2 MeV cm^ g-^, e = 510 GeV, and (3 = 3.92x10-^ cm^ g-^. The cross section 
da/dE** cm^ erg~^ can be expressed by the differential cross section, which is explained 
below, as 

dE^- eJ, dxdy''''- ^^'^ 

The inclusive cross section for the reaction /i^ + N ^ + anything is obtained by 
the renormalization-group-improved parton model (A^ = {n + p)/2 is the isoscalar nucleon). 
In this model, the differential cross section is written in terms of the scaling variables x — 
Q^/2Mv and y^vjE^ as 

2 

{xq{x,Q'')^x{\-vfq{x,Q'')). (63) 

where — is the invariant momentum transfer between the incident neutrino and outgoing 
muon, V = Ey — E^^YS, the energy loss in the laboratory frame, M and M^^ are the nucleon 
and intermediate-boson masses, and is the Fermi constant. As for the distribution 
function of valence and sea quarks, we adopt the fitting formula (Set 2 with A = 290 MeV) 
presented by Eichten-Hinchliffe-Lange-Quigg (EHLQ) for x > 10"^ (Eichten et al. 1984; 
Eichten et al. 1986), while the double- logrithmic approximation (DLA) presented by Reno 
and Quigg (1988) for x < 10^1 



9V _ 2Gj,ME^ 
dxdy TT 



w. 



3. RESULTS 

In this study, when we estimate the flux of neutrinos and gamma-rays from pion decays, 
we have to check whether the energy spectrum of protons can be regarded to obey the rela- 
tivistic Maxwellian distribution. In the cases in which energy loss processes are so effective 
that the energy spectrum of protons do not obey the relativistic Maxwellian distribution any 
longer, the formulation in this study can not be used and Boltzmann equation should be 
needed to estimate the proton's energy distribution and resulting flux of neutrinos, which is 
outscope and next step of this study. Thus, we restrict the situations that meet the following 
constraints in the whole region of the nebula flow. These are: (i) production rate of pions 
[erg s~^] is much smaller than the luminosity of the pulsar wind, (ii) synchrotron cooling 
timescale of protons is longer than traveling timescale (defined below) or timescale of pp 
collisions, (iii) energy transfer timescale from protons to electrons is longer than traveling 
timescale or timescale of pp coUisions (see also the footnotes of table 2) . 
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Here we introduce the definitions of the timescales mentioned above. As for the syn- 
chrotron coohng timescale of protons can be written as 

(in \ 
— ) X te,syn ^ 1-1 X lO^^^e.syn, (64) 

where t^. gyn is the synchrotron energy loss timescale for an electron and can be written as 



. ,9 X 10 (ifX) (151^) (a5) 

Traveling timescale, which means the timescale for the bulk flow of the nebula flow to be 
conveyed to outer region, is defined as 

T 

^travel — -, (66) 
V 

where r is the location of the fluid element measured from the pulsar and v is the bulk speed 
of the nebula flow. Collision timescale between protons is estimated as 

*coi = , (67) 

where n is the number density of protons, app is the cross section of pp interaction, and c is 
the speed of light. For simphcity, to estimate this timescale, app is set to be lOOmb. Collision 
timescale of energy transfer from protons to electrons is estimated as (Stepney 1983) 

4 ne / kT, y 1 

'- = ^-(^1;^ (68) 



In A n„ V TTLpc'^ 



where (Jt is the Thomson cross section and In A is the value of the Coulomb logarithm which 
is estimated as 



In A = In 



hwp 



(69) 



Here ujp is the plasma frequency and is given as 



for a highly relativistic thermal electron gas (Gould 1981). Here we assume that Tg — 
{me/mp)Tp and rie = Up to estimate the energy transfer timescale, although we found that 
the energy transfer timescale depends very weakly on their number ratio. In fact, we found 
that the energy transfer timescale hardly changes even if the number density of electrons 
becomes ~ 10^ times larger than that of protons. 
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Some explanations should be added to the constraints (ii) and (iii). We consider that 
the Maxwellian distribution is hold even if the synchrotron cooling timescale of protons is 
shorter than the pion production timescale, as long as the traveling timescale is shorter than 
the synchrotron cooling timescale. This means the flow is conveyed to outer region without 
losing their energy. It is noted that the definition of the traveling timescale is also frequently 
used as the timescale of adiabatic cooling. However, effect of adiabatic cooling has been 
already taken into consideration because we adopt the formulation of nebula flow. Next, we 
found that the synchrotron cooling timescale of electrons is shorter than the energy transfer 
timescale from protons to electrons whenever the constraint (iii) is not hold. Thus, in the 
cases where the constraint (iii) is not hold, most of the energy of protons is considered to 
be lost through the synchrotron cooling of electrons. Finally, we assume that the energy 
distribution of protons is hold to be Maxwellian. However, exactly speaking, we should 
consider another condition, that is, whether the Coulomb collision timescale between protons 
is sufficiently short that the Maxwellian distribution can be hold. If this condition can not 
be satisfied, temperature of protons can not be defined exactly and formulation of the nebula 
flow must be altered by using Boltzmann equation. However, as explained in section 2.1.2, 
we found in this study that charged and neutral pions are mainly produced just behind the 
termination shock where the energy distribution of protons obey relativistic Maxwellian. 
Thus, we think that the order-estimate of the neutrino flux produced by pp collisions will be 
valid even if the MHD equations are adopted. 

Here we show the results. At flrst, in flgure 5, we show contour of the neutrino event 
whose energy is greater than 10 GeV per year for a km^ detector of high-energy neutrinos 
as functions of age [yr] of a pulsar and initial bulk Lorcnz factor of the pulsar wind (Models 
A2-A9, B2-B9, and C2-C9). The pulsar is assumed to be located 10 kpc away from the 
earth. The amplitude of the magnetic field and period of the pulsar arc assumed to be lO^^G 
and 1ms. It is noted that the spin down age is about 2x10^ yr, so that the range of the age 
is limited to be within 10^ yr. Also, when the age of the pulsar is about 10~^ yr (Models 
D4-D9), synchrotron cooling timescale is so short that the condition (ii) is not satisfied (see 
table 2). As for Model D3, coUision timescale between protons is the fastest one and the 
condition (i) is not satisfied. As for the bulk Lorenz factor of the initial pulsar wind, the 
maximum value is determined by Eq.(28). When the bulk Lorenz factor is about 10^, the 
energy transfer timescale from protons to electrons are too short (see table 2). 

Figure 6 is same with figure 5, but horizontal line represents the initial bulk Lorenz 
factor of pulsar wind with the age of the pulsar fixed as 1 (Models C2-C9), 10^ (Models 
B2-B9), and lO^yr (Models A2-A9). The reason for the event rate to become smaller for the 
case of large initial Lorenz factor is, as shown in figure 7, that the number density of protons 
behind the termination shock becomes smaller. Also, the reason for the event rate to become 
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smaller for the case of the Lorenz factor to be smaller than 10^ is that high energy protons 
do not exist and high energy neutrinos are not produced. It is noted that the cross section 
between neutrino and quark is roughly proportional to energy of neutrinos in the quark-rest 
frame. 

Figure 7 represents the density (sohd line) and temperature (dashed line) behind the 
termination shock as a function of the initial bulk Lorenz factor of the pulsar wind. The 
amplitude of the magnetic field, period, and age of the pulsar are assumed to be lO^^G, 
1ms, and lO^yr (Models A2-A9). As mentioned above, the number density of protons is 
smaller for the larger initial bulk Lorenz factor. We can also confirm that the temperature 
of protons behind the termination shock wave is roughly same with the initial kinetic energy 
of the protons. 

The reason why the flux of neutrinos becomes larger with time can be understood by 
figure 8. Figure 8 shows profiles of velocity, number density of protons, temperature, mag- 
netic field, and emissivity of charged pions in units of [10"^°erg s~^ cm~^] and [10~^°particles 
s~^ cm~^]. The inner and outer boundaries correspond to the location of the termination 
shock wave and the innermost region of the surrounding supernova remnant. The amplitude 
of the magnetic field and period of the pulsar are assumed to be lO^^G and 1ms. The left 
panel represents the case that the age of the pulsar is 1 yr (Model C5), while the right panel 
shows the case that the age of the pulsar is 10^ yr (Model A5). In both cases, initial bulk 
Lorenz factor of the pulsar wind is set to be 10^. The pressure of the innermost region of 
the remnant becomes smaller with time (see figure 2). Thus the location of the termination 
shock moves inward with time, which results in higher density and higher emissivity of neu- 
trinos. Also, additionally, the radius of the remnant becomes outward, which results in the 
region of the nebula flow becomes larger with time. These are the reason why the flux of 
neutrinos becomes larger with time. We can also find that the emissivity of charged pions 
[particles cm~^ s~^] in the nebula fiow can be roughly estimated as 

e = n^ca.. = 3 x 10"^ ( J" , V ( -^J^\ . (71) 

Vl05cm-3y VlOOmb/ ^ ' 

For example, in the left panel of figure 8, the density and emissivity just behind the termi- 
nation shock are about 10^ [cm~^] and 10^'' [particles cm~^ s~^], respectively. Also, those 
values just forward the innermost region of the surrounding supernova remnant are about 
1 [cm~^] and 10~^^ [particles cm~^ s~^], respectively. Thus the emissivity can be roughly 
reproduced by using Eq. (71), although the emissivity is obtained by using formulation pre- 
sented in section 2.2.1. This means the rough estimate holds a good approximation even if 
the protons obey the Maxwellian distribution. Also, as pointed out in section 2.1.2, we can 
find that pions are mainly produced just behind the termination shock. 
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Before we go to the further discussion, we consider the proper value for the bulk Lorenz 
factor of the pulsar wind by introducing the Goldreich- Julian value, which is usually used to 
estimate the charge density of the e^e~ pair plasma, although the initial bulk Lorenz factor 
of the pulsar wind is treated as a free parameter in figure 7. Goldreich- Julian value at the 
equatorial plane at radius r can be written as 

^Gjr = 7^ = 6.9x10-2^ cm-^ 72 

where B{r) is the amplitude of the magnetic field at radius r in a unit of G, O is the angular 
velocity of a pulsar, and P is the period of the pulsar in a unit of sec. Assuming that the 
magnetic field around the pulsar is dipole, Eq.(72) can be written as 

ncj{r) = 6.9 x 10^^ (^-^j (^^J i^—^j cm-^ (73) 

where Bp is the amplitude of the magnetic field at the pole of the pulsar. Since the location 
of the light cyhnder is 

i?ic = 4.8 X lOM I cm, (74) 
\ 1ms / 

the Goldreich- Julian value at the hght cyhnder can be expressed as 

nGa(i?ic) = 6.4 x 10^^ ' cm'^ (75) 

From this value, wc can estimate the net number density of electrons (| n(e+) — n(e~) |) in 
the pulsar wind just ahead of the termination shock as 

f— V^Gj(i?ic) cm-^ (76) 



where we use the conservation law of number flux of protons, AnRfjI'nGjiRic) — 47rr^rn(rs). 
Eq.(76) can be written by inserting Eq.(75) as 

ImsV f Bp \ fW^ cm 



n(r.) = 1.6 X 10 ^— j j ^— j cm-. (77) 

The location of the termination shock is found to be ~ 10^^ cm from figure 8. On the other 
hand, the number density of protons can be estimated by using relation L ~ A^nnV'^r^mpC^ 
(see Eqs.(2) and (29). note that a is much smaller than unity and bulk flow is highly 
relativistic) . It can be expressed as 
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By comparing Eq.(77) with Eq.(78) and assuming that number density of proton is compa- 
rable with the net number density of electrons (| n{e'^) — n{e~) |), proper value for the bulk 
Lorenz factor seems to be F ~ 10^ - 10^. 

In figure 9, we show spectrum of energy fluxes of neutrinos from a pulsar which is located 
lOkpc away from the earth. The amphtude of the magnetic field and period of the pulsar 

is assumed to be lO^^G and 1ms. The detection limits of the energy flux for AMANDA- 
BIO, AMANDA II (lyr), and IceCube are represented by horizontal lines. The atmospheric 
neutrino energy fluxes for a circular patch of 1° are also shown (Honda et al. 1995). The 
left panel represents the case that the age of the pulsar is lyr (Models C2-C9), while right 
panel represent the case that the age is lO^yr (Models A2-A9). From this figure, we can find 
that there is a possibility to detect the signals of neutrinos from pulsar winds in our galaxy. 

In figure 10, we show neutrino event rate per year from the nebula fiow considered in 
figure 9 as a function of the muon energy threshold in a km^ high-energy neutrino detector. 
The left panel represents the case that the age of the pulsar is lyr (Models C2-C9), while 
right panel represent the case that the age is lO^yr (Models A2-A9). The event rates expected 
of atmospheric neutrino for a circular patch of 1° are also shown (Honda et al. 1995). The 
event rates are shown for the vertical downgoing atmospheric neutrinos and horizontal ones. 
It is confirmed that there is a possibility that the signals from the nebula flow dominates 
the background of atmospheric neutrinos. 

For comparison, in flgure 11, we show the spectrum of energy fluxes of neutrinos from 
a pulsar of which the amplitude of the magnetic held and period are assumed to be lO^^G 
and 5ms. Location of the pulsar is assumed to be 10 kpc away from the earth. The left 
panel represents the case that the age of the pulsar is lOyr (Models G1-G8), while right 
panel represents the case that the age is lO^yr (Models E1-E8). In these cases, the fluxes 
of neutrinos are much lower than the atmospheric neutrinos and detection limits of km^ 
high-energy neutrino detectors. We can conclude that the detectability of the signals from 
the pulsar winds strongly depends on the period of the pulsar, which reflects the luminosity 
of the pulsar winds (see Eq.(29)). 

We check the reason why the flux of neutrinos becomes lower for the case where period 
of the pulsar is 5ms. Figure 12 shows proflles of velocity, number density of protons, tem- 
perature, magnetic fleld, and emissivity of charged pions in units of [lO^^'^erg s~^ cm~^] and 
[10~^°particles s~^ cm~^]. The left panel represents the case that the period and age of the 
pulsar are 1ms and 10^ yr (Model A5; same with the right panel of flgure 8), while the right 
panel shows the case that the period and age of the pulsar are 5ms and 10^ yr (Model E5). 
In both cases, initial bulk Lorenz factor of the pulsar wind is set to be 10^. In the latter 
case, the location of the termination shock must be outward compared with the former in 
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order to achieve the pressure balance between the outermost region of the nebula flow and 
innermost region of the remnant. This is because luminosity of the pulsar wind and the 
pressure behind the termination shock are lower in the latter case. That is why the density 
of behind the termination shock is lower in the latter case and results in lower emissivity of 
neutrinos. 

Here we consider the detectability of gamma-rays from neutral pions. In figure 13, 
integrated gamma-ray fluxes from the neutral pion decays are shown assuming that the 
supernova ejccta has been optically thin for gamma-rays (Matz et al 1988). The amplitude 
of the magnetic field and period of the pulsar arc assumed to be lO^^G and 1ms. The left 
panel represents the case that the age of the pulsar is lyr (Models C2-C9), while the right 
panel shows the case that the age of the pulsar is 10^ yr (Models A2-A9). The detection limits 
of integrated fluxes for GLAST, STACEE, CELESTE, HEGRA, CANGAROO, MAGIC, 
VERITAS, and H.E.S.S. are also shown. Prom these figures, we can find that there is a 
possibility to detect gamma-rays from decays of neutral pions by these telescopes. As for 
the spectrum of gamma-rays from electrons, it is difficult to estimate it because we have 
to solve time evolution of energy distribution of electrons taking into account the effects of 
cooling processes such as synchrotron emission, inverse compton, and pair annihilation as 
well as heating process such as Coulomb collisions between electrons and protons and pair 
creations. This is difficult to solve and we regard it as a future work. 

On the other hand, we show the same figures in figure 14 but for the case that the 
period of the pulsar is 5ms. The left panel represents the case that the age of the pulsar is 
lOyr (Models G1-G8), while the right panel shows the case that the age of the pulsar is 10^ 
yr (Models E1-E8). In these cases, the flux of gamma-rays is too low to detect. 

Finally, it is noted that there are some cases in which the location of the termination 
shock is driven back to the surface of the neutron star since the pressure balance can not be 
achieved (see section 2.1.3). Models H2 and H3 are the case (see table 2). In these cases, the 
density behind the shock is so large that the resulting flux of neutrino is very high. There 
will be a tendency that the reverse shock is hkely to be driven back to the surface of the 
neutron star when the age of the pulsar is young, because the pressure at the innermost 
region of the remnant is large (see figure 2). Thus there may be a possibility to detect a 
number of high energy neutrinos from the pulsar wind at the early stage of the supernova 
explosion, although the protons in the nebula flow will also suffer from synchrotron cooling 
and/or energy transfer to electrons. 
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4. DISCUSSIONS 

First, we discuss the timescale of energy loss of protons due to the inverse compton (IC) 
scattering. This timescale depends on the strength of the seed photons that suffers from 
great uncertainty. However, we can estimate the lower limit of this timescale and show that 
this process is not so effective. The ratio of the synchrotron coohng timescale of protons 
relative to IC can be written as 

^sync _ Uj 

where and Ub are energy density of seed photon and electric- magnetic fields. In the 
case in which energy transfer from protons to electrons is ineffective, the energy density of 
seed photon should be smaller than that of electrons as long as the seed photons come from 
bremsstrahlung and/or synchrotron emission from electrons. Thus, using Eq.(22), the upper 
limit of the ratio of tgync/^ic can be expressed as 

he ~ Ub~ UpTUpUB ~ npTUp a{3az^ + (3a) 2/3^4/3 14/3 [4 + + {3a/ zy/^}^] ' ^ ^ 
We show in figure 15 the ratio of the thermal energy of electrons relative to the energy of 
electric- magnetic fields as a function of z in the case where Ue—rip. Here a is set to be 
(Tc = 6.67 X 10~^. We can find that the upper limit of the ratio of igync/^ic is smaller than 
unity as long as Ue/up < 2 x 10^. That is why we can conclude that IC cooling is not so 
effective. Here we neglect the contribution of cosmic microwave-background, since it can be 
expressed by an equivalent field strength i?cMB = 3.24/iG, which is sufficiently small and can 
be neglected. We also check the contribution of the photon field in the supernova remnant 
around the nebula fiow, since a portion of it may be transported into the region of the nebula 
fiow. The energy density of the photon field in the supernova remnant can be estimated from 
figure 2 and expressed by equivalent field strength 5sn ~ 3.9 G, 1.3x10"^ G, 4.1x10"^ G, 
and 1.6x10"^ G ior t — 1, 10, 10^, and 10^ yr. These values are also smaller than strength 
of the magnetic fields in almost all cases (see figures 8 and 12). We should also estimate 
the number density of photons from surface of the pulsar. The temperature of the surface 
of the neutrons star is considered to become in the range (10^ - 10^) K about 1 yr after the 
explosion, although there are some uncertainties of the cooling processes (Slane et al. 2002). 
Thus the energy density of the photon field at radius r can be estimated as 

_47rRVT4 

(^.pulsar — ^^^2^ ' vol) 

where R is the radius of the pulsar and a is the Stefan-Boltzmann constant. Eq. (81) can 

be expressed by the equivalent field strength as 

/ R \ r T \ /IO12 m\ 
^Pu.. = 2.2xl0(^j(^j(^j G, (82) 
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which is also smaller than strength of the magnetic fields in almost all cases (see figures 8 
and 12). 

Next, we consider the effect of synchrotron cooling of muons and pions. As for the mean 
fife time of muon in the rest frame is 2.197x10"^ s and that of charged pion is 2.603x10"^ 
s. Since synchrotron cooling time of muon is 



the condition that the synchrotron cooling time becomes shorter than that of the mean life 
time in the fiuid rest frame can be expressed as 



Since the amplitude of the magnetic field in the nebula flow becomes ~ 10^ G at most (see 
figures 8 and 12), the energy spectrum of neutrinos whose energy is larger than ~ 10^ GcV 
may be modified due to the efi^ect of the synchrotron cooling of muons. Similarly, the energy 
spectrum of neutrinos whose energy is larger than ~ lO'' GeV may be modified by the effect 
of synchrotron cooling of charged pions. 

Third, we discuss the cooling process due to photopion production which is interesting 
and should be investigated in detail in the forth-coming paper because this is another process 
of producing pions and may make the flux of neutrinos and gamma-rays enhance. Here we 
give the rough estimation for this effect. The source of the seed photons are considered to 
be the synchrotron photons emitted by electrons, the photon fleld in the supernova remnant 
around the nebula flow, and photons from surface of the pulsar. The typical frequency of 
the synchrotron photon is (Rybicki and Lightman 1979) 



where 7 is the Lorenz factor of electrons in the fluid rest frame and is considered to be almost 
same with the initial bulk Lorenz factor of the pulsar wind (= F) as long as the electrons do 
not lose their energy by synchrotron cooling. Here the average angle between the magnetic 
fleld and direction of motion of electron is set to be 7r/4. This frequency corresponds to the 
energy 



t. 




(83) 




(84) 




(85) 
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Since there is a large peak in the photopion production cross section at e ~ 0.35mpC^ at the 
rest frame of a proton and the width of the peak cross section is also Ae ~ 0.2GeV (Stecker 1979), 
the energy of protons in the fluid rest frame which interact with the seed photons can be 
expressed as 

0.35m^c^ „ /10^G\ , , 

E, ^ = 1.2 X 10^^-2 (^-^J GeV. (87) 

This shows, to be sure, that the protons in the nebula flow cause the photopion interactions 
when the initial bulk Lorenz factor of the pulsar wind is larger than ~ 10^. We roughly esti- 
mate the upper limit of the contribution of photopion production relative to pp interaction. 
The maximum number density of photons by synchrotron emission is estimated as 

^Ue_ _ neirie Up 

7,sync _ ^^^^ 

~ — - = 2.0 X loM — — — He. 88 

Thus, if we assume that the energy loss rate per interaction is same between photopion 
production and pp interaction, and assume that the cross section of photopion production is 
constant (~ 5 x 10~^^cm^) , the contribution of photopion production relative to pp interaction 
(~ 10~^^cm^) can be estimated to be 

^ 1.0 X 10^ f V!^V (89) 

In reality, as mentioned above, there is a large peak in the photopion production cross 
section at the photon energy e = 0.35mpC^ in the proton rest frame due to the A(1232) 
resonance (Stcckcr 1979; Waxman and Bahcall 1997; Amato et al. 2003). Thus the contri- 
bution of photopion production would be much smaller than the estimated value in Eq.(89) 
in the case that the typical energy of protons which interact with the synchrotron photons 
from electrons through the A(1232) resonance is much different from the mean energy of 
protons in the fluid. We emphasize again that we estimated here just the upper limit of the 
contribution of photopion production. Also, it is noted here that the synchrotron cooling 
time of electrons can be estimated as 

te,sy„c = 7.6 X 10-^ (H!Ul^y 3. (90) 



Thus, even if the effect of the photopion production can not be neglected, it will be effective 
only within a very short distance from the location of the termination shock. As for the 
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contribution of the photon field from the supernova remnant, the energy of protons that 
interacts with the photon field can be estimated as 



where Ty is the temperature of the photon field. Thus the required energy for a proton to 
interact with the photon field from the supernova remnant seems to be too high. As for the 
contribution of the photon field from the pulsar, the number density of the photon field can 
be expressed as 



Thus, from Eqs.(88) and (92), the photon fields from the pulsar can dominate the synchrotron 
photons. As for the energy of protons that interacts with the photon field, it can be estimated 

as 



Thus, protons with energies ~ 10^ GeV may interact with the photons from the pulsar and 
produce much more pions than estimated in this study. Thus, the spectrum of neutrinos 
and gamma-rays around ~ 10^ GeV may be modified and intensity may be enhanced due 
to this effect. This effect is very important and should be investigated carefully in the very 
near future. 

In this study, we consider the pion production due to pp interaction. However, in many 
cases, total energy of colliding protons in their center-of-mass system exceeds the experi- 
mental limit. Moreover, it exceeds the mass of the top quark {rrit = 174.3 ± S.lGeV from 
direct observation of top events, and rrit = ITS.llg^g'^GeV from standard model electroweak 
fit (Hagiwara et al. 2002)). As explained in section 2.2.2, it is confirmed that the scahng 
assumption does not much affect the resulting gamma-ray spectrum by comparing the result- 
ing gamma-ray spectrum up to ~ 10^ GeV with the gamma-ray production by the mini-jet 
model including the QGD effect even if the scaling model is extrapolate to higher energies. 
However, it should be confirmed whether the scaling assumption is still valid at more high 
energies, which is investigated in the near future. 

In this study, we have estimated the energy transfer timescale assuming that it is mainly 
caused by Goulomb scattering between protons and electrons (Stepney 1983). In this case. 
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the energy transfer timescale depends very weakly on the number ratio of protons to electrons 
(positrons), as shown in section 3. However, it is pointed out by Hoshino et al. (1992) that 
the efficiency of the energy transfer from protons to positrons may be higher than that 
from protons to electrons due to the effect of cyclotron resonance. Thus, if this process is 
really effective and determines the energy transfer rate, number ratio of protons to electrons 
(positrons) should be important because the energy transfer timescale will strongly depend on 
the number ratio. Also, the expected flux of neutrinos and gamma-rays from pp interactions 
will be decreased than estimated in this study. In this case, however, emission of gamma- 
rays from electron-positron plasma should be important through the process of synchrotron, 
inverse compton, and pair annihilation. Neutrino emission from electron-positron plasma 
is also expected from electron-positron pair annihilation. In more realistic cases, there will 
be heavy nuclei in the pulsar wind, which should be taken into account as a next step of 
this study. In this case, hadronic interactions should be enhanced since heavy nuclei are 
composed of protons and neutrons and baryon number should be enhanced. In order to 
obtain the contribution of heavy nuclei, photodisintegration of heavy nuclei should be also 
taken into account, which will make the situation rather complicated. 

In this study, we assumed that the pressure balance is hold at the interface between 
the nebula flow and supernova remnant. Also, we fixed the fraction of magnetic field energy 
in the pulsar wind so that the velocity of the nebula flow is almost same with that of the 
supernova remnant at the interface, which makes the location of the termination shock fixed. 
However, in reality, the pressure balance should not be always hold and velocity of the nebula 
flow and supernova remnant will be different in general. This means that we need to consider 
the momentum transfer between the nebula flow and supernova remnant, which will change 
profiles of density, temperature, and amplitude of the magnetic field in the nebula flow. 
This is a future work and it may require numerical calculations. Also, it is assumed that 
the system is spherical symmetry for simplicity. Of course, in reality, we should consider the 
effects of geometrical anisotropy of the system (Hester et al. 2002). 

In this study, we investigated the effects that bulk Lorenz factor is dissipated into the 
random motion of the charged particles when they pass through the termination shock. 
Of course, a part of charged particles will be accelerated at more high energies due to 
Fermi acceleration (Blandford and Ostriker 1978). Maximum energy of protons due to Fermi 
acceleration can be estimated using Larmor radius as 

B_ = 3.3xlO»(j|^)(^)GeV. (94) 

From flgures 8 and 12, E^ax will be about 10^ GeV and 10^ GeV for the cases of P = 1ms 
and 5ms {Bp is set to be lO^^G). Thus, the flux and mean energies of neutrinos and gamma- 
rays may be enhanced due to the effects of Fermi acceleration. We will have to investigate 



31 



carefully how much of charged particles are accelerated by Fermi acceleration in order to 
estimate this effects quantitatively. 

It is noted that the is generated due to the effect of neutrino oscillation (Fukuda et al. 1998). 
The scale length for with energy Ejj to be converted to be v,,- due to the vacuum oscillation 
is written as 



where Am is the mass difference between m2 and 777.3. Thus the neutrinos produced in the 
nebula flow are well mixed and ratios of the flux among three flavors (e, //., and r) will 
become order unity at the earth. Thus there may be a possibility to detect the double-bang 
events (a big hadronic shower from the initial Vt interaction and a second big particle cascade 
due to T decay, which was pointed out by Learned and Pakvasa 1995) at the km^ neutrino 
telescopes from the system investigated in this study. 

At present, there is no candidate of a pulsar whose wind is more luminous than ~ 10^^ 
erg s"^ like the case ol B — lO^^G and P — 1ms in our galaxy (Manchester et al. 1999; 
Torres and Nuza 2003). This is consistent with the fact that there is no point- like source 
detected at AMANDA-BIO detector (Ahrens ct al. 2003). We hope that such point-hke 
sources are found at AMANDA-II and/or IccCubc detectors in the near future. Also, we 
hope that such an active pulsar is found by pulsar surveys in radio. X-ray, and gamma-ray 
bands. We also hope that a newly born supernova will appear in our galaxy in the near 
future. 

Recently, Granot and Guetta (2003) pointed out the possibility to detect neutrinos 
from pulsar nebula that is associated with a gamma-ray burst (GRB). It will be interest- 
ing to apply our model to GRB scenario and estimate fluxes of neutrino background by 
setting the amplitude of magnetic held at the pole of the pulsar is set to be (10^^ - 10^^) 
G (Granot and Guetta 2003). Information on GRB formation rate history and explosion 
mechanism of GRB may be obtained if such signals are detected (Nagataki et al. 2003a; 
Nagataki et al. 2003b). 



In this study, we have estimated fluxes of neutrinos and gamma-rays that are generated 
from decays of charged and neutral pions from a pulsar surrounded by supernova ejecta in our 
galaxy, including an effect that has not been taken into consideration, that is, interactions 
between high energy cosmic rays themselves in the nebula flow, assuming that hadronic 
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SUMMARY AND CONCLUSION 
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components be the energetically dominant species in the pulsar wind. Bulk flow is assumed 
to be randomized by passing through the termination shock and energy distribution functions 
of protons and electrons behind the termination shock obey the relativistic Maxwellians. 

We have found that fluxes of neutrinos and gamma-rays depend very sensitively on the 
wind luminosity, which is assumed to be comparable with the spin-down luminosity. In 

the case where B = lO^^G and P = 1ms, neutrinos should be detected by km^ high-energy 
neutrino detectors such as AMANDA and IceCubc. Also, gamma-rays should be detected by 
Cherenkov telescopes such as CANGAROO and H.E.S.S. as well as by gamma-ray satellites 
such as GLAST. On the other hand, in the case where B = lO^^G and P = 5ms, fluxes of 
neutrinos will be too low to be detected even by the next-generation detectors. However, 
even in the case where B — lO^^G and P — 5ms, there is a possibihty that very high flux 
of neutrinos and gamma-rays may be realized at the early stage of a supernova explosion 
{t < lyr) , where the location of the termination shock is very near to the pulsar (ex. Model 
H3 and H4 in table 2). We also found that there is a possibility that protons with energies 
~ 10^ GeV in the nebula flow may interact with the photon field from surface of the pulsar 
and produce much pions, which enhances the intensity of resulting neutrinos and gamma- 
rays. 

We have found that interactions between high energy cosmic rays themselves are so 
effective that this effect can be confirmed by future observations. Thus, we conclude that it 
is worth while investigating this effect further in the near future. 

The author thanks K. Sato, M. Hoshino, S. Shibata, S. Yamada, S. Yoshida, T. Naito, 
and K. Kohri for useful comments. The author is also pleased to acknowledge useful com- 
ments by the anonymous referee of this manuscript. This work is supported in part by 
Grants-in-Aid for scientific research provided by the Ministry of Education, Science and 
Culture of Japan through Research Grant No.S 14102004 and No. 14079202. 



A. Estimation of the Adiabatic Index 



According to Hoshino et al. (1992), the distribution functions of protons behind the 
termination shock are almost exactly described by relativistic Maxwellians as 



iV(7) = A'jexp 



rripC 



(7-1) 



(Al) 



where is the Boltzmann constant and A is the normalization factor in units of cm ^. In 
this case, the number density (n) [cm~^], energy density (e) [erg cm~^], and pressure (P) 
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[dyn cm ^] can be written as 



1 1 ■ 
- + — 



(A2) 



2 ^pC7 



Nil) , 



'2 2 
1 + - + — 



(A3) 



and 



1 /■ / rfn \ ,o I [ pc^ iV(7) 2 , Am„c2 „ . , -HiEf!^ , 
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1 1 ■ 
- + — 



(A4) 



where a — nipC^ /kBT, E — nipC^i, dn/d^xd^p is the number density in phase space. In this 
case, the relation P — 2nkBT/3 holds exactly, and in the high temperature hmit (a <^ 1), 
the relation 



E = eV = 3PV 



(A5) 



holds approximately, where V is the volume [cm^] . In this case, the relation 3VdP — —APdV 
holds in the adiabatic transition. Thus P oc V~^^^ oc p'^/^ and adiabatic index becomes 4/3 
in high temperature limit. 



B. Derivation of Number Flux of Pions in Observer's Frame 



We consider an fluid element and number spectrum of pions produced in the fluid 
element. Using the number spectrum of pions per solid angle [particles s~^ erg~^ sr~^], the 
conservation law of number of particles can be expressed as 

dE^ f sin eded^-^F{E^)dtT{l -pcosO), (Bl) 



dE^ / sin eded(p—F{E^)dt 



where bars are labeled for the quantum of the observer's frame, while the quantum without 
label denote the ones in the fluid-rest frame. The axis ^ = ^ = is set to be ahgned with 
the direction of the observer (i.e. the earth) measured from the fluid element, f is the bulk 
Lorenz factor of the fluid element in the observer's frame. It is noted that dtA — r{l — P)dt 
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is the time interval of the radiation as received by a stationary receiver in the observer's 
frame (Rybicki and Lightman 1979). Here we define the quantum N{9) as 



J dE^F{E^) [particles s"^]. (B2) 



Using this expression, the relation 

dQ r(l-/5cos^)sin^#47r ^ ^ 

can be derived. Using the relation cos^ = (cos 6* + + PcosO) and d4> = d4>, Eq. (B3) 

can be expressed as 

Also, the relation F{E^) — F{E^)dEi,/dEi, can be expressed as 

Finally, using Eqs. (B4) and (B5), the relation 

-Lf(E,) = i = — -F(E,) (B6) 

dQ ^ r2(l-/3cos^^)2 47r ^ ^ ' 

is derived. 



C. Analytical Solution for the Nebula Flow 



The analytical solution of equations for the nebula flow is obtained by Kennel and 
Coroniti (1984), which is expressed as 



;i + «y)V2 



5 + A(^;;22)-i/3 + i 



r2(l + 5 + A), 



(CI) 



where v is defined d&v — ujui^ z\% defined as 2; = r/r^, 5 is defined as 5 = 47rn2r|mpC^/ B\ ~ 
u^juxu^ and A is defined as \^'kP'^\IB\. The solution for this equation can be written as 

' (C2) 



2(1 + A) 



where 



G = 1 + {1 + + [(1 + x'^f - 1] + {1 + - [(1 + x-^f - 1] 



(C3) 
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X 



z/za, 



(C4) 



and 



za 



'2 ' 
_27(1 + A)2_ 



1/2 



(C5) 



Using this solution, we can derive the solution for the total pressure Pt in the postshock 
in Eq. (22). 



REFERENCES 

Ahrens, J., Bai, X., Barouch, G., Barwick, S. W., Bay, R. C, Becka, T., Becker, K.-H., 
Bertrand, D., Binon, F., et al., 2003, ApJ, 583, 1040 

Amato, E., Guetta, D., Blasi, P., 2003 A&A, 402, 827 

Atoyan, A. M., Aharonian, F. A., 1996, MNRAS, 278, 525, (1996) 

Badhwar, G. D., Stephens, S. A., Golden, R. L., 1977, Phys. Rev. D, 15, 820 

Beall, J. H., Bednarek, W., 2002, ApJ, 569, 343 

Bednarek, W., Prhotheroe, R. J., 1997, Phys. Rev. Lett., 79, 2616 

Bednarek, W., 2001, A&A, 378, 49 

Bednarek, W., Bartosik, M, 2003, A&A, 405, 689 

Bednarek, W., 2003, astro-ph/0305430 

Bednarek, W., 2003, astro-ph/0307216 

Berezinsky, V. S., Gaisser, T. K., Halzen, F., Stanev, T., 1993, Astropart. Phys., 1, 281 
Blasi, P., Epstein, R. I., Olinto, A. V., 2000, ApJ, 533, 123 
Blandford, R. D., Ostriker, J. P., 1978, ApJ, 221, L29 
Dermer, C. D., 1986a, ApJ, 307, 47 
Dermer, C. D., 1986b, A&A, 157, 223 

Eichten, E., Hinchliffe, I., Lane, K., Quigg, C., 1984, Rev. Mod. Phys., 56, 579 



-36- 



Eichten, E., Hinchliffe, I., Lane, K., Quigg, C, 1986, Rev. Mod. Phys., 58, 1065 

Fukuda, Y. et al. (Super-Kamiokande Collaboration) Phys. Rev. Lett. 81, 1562 (1998) 

Gaisser, T. K., Grillo, A. P., 1987, Phys. Rev. D, 36, 2752 

Gaisser, T. K., Halzen, P., 1987, Phys. Rev. Lett., 54, 1754 

Granot, J., Guetta, D., 2003, Phys. Rev. Lett., 19, 1102 

Goldreich, P., Juhan, W. H., 1969, ApJ, 157, 869 

Gould, R. J., 1981, Phys. Fluids 24, 102 

Gunn, J. E., Ostriker, J. P., 1969, Phys. Rev. Lett., 22, 728 

Haas, R., Colgan, J., Erickson, P., Lord, D., Burton, G., HoUenbach, J., 1990, ApJ, 360, 257 

Hagiwara, K. Hikasa, K., Nakamura, K., Tanabashi, M., Aguilar-Benitez, M., Amsler, C., 
Barnctt, R. M., Burchat, P. R., Carone, C. D., Caso, et al., 2002, Phys. Rev. D, 66, 
010001 

Hashimoto, M., 1995, Prog. Theor. Phys., 94, 663 

Hester, J. J., Mori, K., Burrows, D., Gallagher, J. S., Graham, J. R., Halverson, M., Kader, 
A., Michel, F. C., Scowen, P., 2002, ApJ, 577, L49 

Honda, M., Kajita, T., Kasahara, K., Midorikawa, S., 1995, Phys. Rev. D, 52, 4985 

Hoshino, M., Arons, J., Gallant, Y. A., Langdon, A. B., 1992, ApJ, 390, 454 

Kennel C. P., Coroniti, F. V., 1984, ApJ, 283, 694 

Landau, L. D., Lifshitz, E. M., 1975, The Classical Theory of Fields, 4th Rd., Pergamon, 
Oxford 

Learned, J. G., Pakvasa, S., 1995, Astropart. Phys., 3, 267 
Mahadevan, R., Narayan, R., Kronik, J., 1997, ApJ, 486, 268 

Manchester, R. N., Lync, A. G., Camilo, P., Kaspi, V. M., Stairs, 1. H., Crawford, P., Morris, 
D., J., Bell, J. P., Amico, N. D., 1999, astro-ph/9911319 

Matz, S. M., Share, G. H., Leising, M. D., Chupp, E. L., Vestrand, W. T., Purcell, W. R., 
Strickman, M. S., Reppin, C, 1988, Nature, 331, 416 



-37- 

Nagataki, S., Hashimoto, M., Sato, K., Yamada, S., 1997, ApJ, 486, 1026 
Nagataki, S., 2000, ApJS, 127, 141 

Nagataki, S., Kohri, K., Ando, S., Sato, K., 2003, Astropart. Phys., 18, 551 

Nagataki, S., Mizuta, A., Yamada, S., Takabe, H., Sato, K., 2003, ApJ, in press. 

Naito, T., Takahara, F., 1994, J. Phys. G: Nucl. Part. Phys., 20, 477 

Reno, M. H., Quigg, C, 1988, Phys. Rev. D, 37, 657 

Ruderman, M. A., Sutherland, P. G., 1975, ApJ, 196, 51 

Protheroe, R. J., Bednarek, W., Luo, Q., 1998, Astropart. Phys., 9, 1 

Rees, M. J., Gunn, J. E., 1974, MNRAS, 167, 1 

Rybicki, G. B., Lightman, A. P., 1979, Radiative Process in Astrophysics (New York: John 
Wiley & Sons) 

Shibata, S., 1991, ApJ, 378, 239 

Slane, P., Helfand D. J., Murray, S., 2002, ApJ, 571, L45 
Spyromilio, J., Meikle, S., Allen, A., 1990, MNRAS, 242, 669 
Stecker, F. W., 1971, Cosmic Gamma Rays, NASA SP-249 
Stecker, F. W., 1979, ApJ, 228, 919 
Stephens, S. A., Badhwar, G. D., 1981, Ap&SS, 76, 213 
Stepney, S., 1983, MNRAS, 202, 467 

Thielemann, F.-K., Nomoto, K., Hashimoto, M., 1996, 460, 408 
Torres, D. F., Nuza, S. E., 2003, ApJ, 583, L25 
Waxman, E., Bahcall, J. N., 1997, Phys. Rev. Lett., 78, 2292 
Woosley, S. E., Weaver, T. A., 1995, ApJS, 101, 181 



This preprint was prepared with the AAS macros v5.0. 



-38- 




Fig. 1. — Solid line: spin down age [yr] as a function of period of a pulsar. Spin down age is 
defined as the time interval for a pulsar's angular velocity to become half of the initial one. 
Dashed line: spin down luminosity as a function of a period of a pulsar. In this study, wind 
luminosity is assumed to be comparable with the spin down luminosity. Amplitude of the 
magnetic field at the pole of a pulsar is set to be lO^^G throughout in this study. 
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Fig. 2. — Relation between pressure [dyn cm~^] and age [yr] of the innermost region of the 
remnant. The discontinuity at t ~ 5 yr reflects the transition from photon-dominated phase 
to matter-dominated phase. This happens when the optical depth of the supernova ejecta 
becomes lower than unity and pressure of photon fields is set to be zero. 
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Fig. 3. — Total pressure as a function oi z = r/rg, normalized by 2ni7iMimpC^/3. a is set to 
be 6.7 X 10-=^ = (7c. 
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Fig. 4. — Sketch of the geometry concerning the individual scattering events. The left panel 
shows the geometry in the fluid-rest frame, while the right panel shows the geometry in the 
rest frame of particle 1. 
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Fig. 5. — Contour of the neutrino event whose energy is greater than 10 GeV per year for 
a km^ detector of high-energy neutrinos as functions of age [yr] of a pulsar and initial bulk 
Lorenz factor of the pulsar wind. The pulsar is assumed to be located 10 kpc away from 
the earth. The amplitude of the magnetic field and period of the pulsar are assumed to be 
10^^ G and 1ms. It is noted that the spin down age is about 2x10^ yr, so that the range of 
the age is limited to be within 10^ yr. Also, when the age of the pulsar is about 10~^ yr, 
synchrotron coohng timescale and/or energy transfer timcscalc from protons to electrons are 
so short that the Maxwcllian distribution of protons is distorted and production of pions is 
severely limited (see table 2). As for the bulk Lorenz factor of the initial pulsar wind, the 
maximum value is determined by Eq.(28). When the bulk Lorenz factor is about 10^, the 
energy transfer timescale from protons to electrons are too short (see table 2) . 
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Fig. 6. — Same with figure 5, but horizontal line represents the initial bulk Lorenz factor 
of pulsar wind with the age of the pulsar fixed as 1 (Models C2-C9), 10^ (Models B2-B9), 
and lO^yr (Models A2-A9). The reason for the event rate to become smaller for the case 
of large initial Lorenz factor is that the number density of protons behind the termination 
shock wave becomes smaller (see figure 7). Also, the reason for the event rate to become 
smaller for the case of the Lorenz factor to be smaller than 10^ is that high energy protons 
do not exist and high energy neutrinos are not produced. It is noted that the cross section 
between neutrino and quark is roughly proportional to energy of neutrinos in the quark-rest 
frame. 
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Fig. 7. — The density (solid line) and temperature (dashed line) behind the termination 
shock wave as a function of the initial bulk Lorenz factor of the pulsar wind. The amphtude 
of the magnetic field, period, and age of the pulsar are assumed to be lO^^G, 1ms, and lO^yr 
(Models A2-A9). 
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Fig. 8. — Profiles of velocity, number density of protons, temperature, magnetic field, and 
emissivity of charged pions in units of [10~^°erg s~^ cm~^] and [10~^°particles s~^ cm~^]. The 
inner and outer boundaries correspond to the location of the termination shock wave and 
the innermost region of the surrounding supernova remnant. The amplitude of the magnetic 
field and period of the pulsar are assumed to be lO^^G and 1ms. Initial bulk Lorenz factor 
of the pulsar wind is set to be 10^. The left panel represents the case that the age of the 
pulsar is 1 yr (Models C5), while the right panel shows the case that the age of the pulsar 
is 10^ yr (Models A5). 
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Fig. 9. — Spectrum of energy fluxes of neutrinos from a pulsar which is located 10 kpc away 
from the earth. The amplitude of the magnetic field and period of the pulsar is assumed 
to be lO^^G and 1ms. The minimum detectable energy flux of AMANDA-BIO, AMANDA 
11 (lyr), and IceCube is represented by horizontal lines. The atmospheric neutrino energy 
fluxes for a circular patch of 1° arc also shown. The left panel represents the case that the 
age of the pulsar is lyr (Models C2-C9), while right panel represent the case that the age is 
lOVr (Models A2-A9). 
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Fig. 10. — Neutrino event rate per year from a pulsar wind as a function of the muon energy 
threshold in a km^ high-energy neutrino detector. The pulsar is assumed to be located 10 
kpc away from the earth. The amplitude of the magnetic field and period of the pulsar is 
assumed to be lO^^G and 1ms. The left panel represents the case that the age of the pulsar 
is lyr (Models C2-C9), while right panel represents the case that the age is lOVr (Models 
A2-A9). The event rates expected of atmospheric neutrino for a circular patch of 1° are also 
shown. The event rates are shown for the vertical downgoing atmospheric neutrinos and 
horizontal ones. 
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Fig. 11. — Same as figure 9, but for the period of the pulsar is 5ms. The left panel represents 
the case that the age of the pulsar is lOyr (Models G1-G8), while right panel represent the 
case that the age is lO^yr (Models E1-E8). Note that the scale of the vertical line is different 
from that of figure 9. 
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Fig. 12. — Profiles of velocity, number density of protons, temperature, magnetic field, and 
emissivity of charged pions in units of [10^^'' erg s^^ cm^^] and [lO^^^'particles s^^ cm^'^]. The 
inner and outer boundaries correspond to the location of the termination shock wave and 
the innermost region of the surrounding supernova remnant. The amplitude of the magnetic 
field is assumed to be lO^^G. Initial bulk Lorenz factor of the pulsar wind is set to be 10^. 
The left panel represents the case that the period and age of the pulsar are 1ms and 10^ yr 
(Model A5; same with the right panel of figure 8), while the right panel shows the case that 
the period and age of the pulsar are 5ms and 10^ yr (Model E5). 
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Fig. 13. — Integrated gamma-ray fiuxes from the neutral pion decays are shown. The 
amplitude of the magnetic field and period of the pulsar are assumed to be lO^^G and 
1ms. The left panel represents the case that the age of the pulsar is lyr (Models C2- 
C9), while the right panel shows the case that the age of the pulsar is 10^ yr (Models A2- 
A9). The minimum detectable integrated fluxes of GLAST, STAGEE, GELESTE, HEGRA, 
GANGAROO, MAGIG, VERITAS, and H.E.S.S. are also shown. 
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Fig. 14. — Same with figure 13, but for the case that the period of the pulsar is 5ms. The 
left panel represents the case that the age of the pulsar is lOyr (Models G1-G8), while the 
right panel shows the case that the age of the pulsar is 10^ yr (Models E1-E8). 
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Fig. 15. — Ratio of the thermal energy of electrons relative to the energy of electric- magnetic 
fields as a function of z (= r/r^) in the case where Ue—rip. Definition of a is shown in Eq.(l) 
and set to be Uc — 6.67 x 10~^. 



-50- 



Table 1. Parameters for the representation of invariant cross section 



Particle 


A 


B 


r 




C 


C 




[mb/(GeVVc3)] 


[(GeV/c)-^] 






[(GeV/c)-^] 


[(GeV/c)-^] 


tt" 
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5.43 


2 


6.1 


-3.3 


0.6 


7r+ 


153 


5.55 


1 


5.37 


-3.5 


0.83 


7r~ 


127 


5.30 


3 


7.03 


-4.5 


1.67 



-51 - 



Table 2. Input and Output Parameters 



Model 
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Model 


E5 


IQi-''' 


5 


1.5 


X 10*0 


10^^ 


10"' 


2.8 X 


10" 


^travel 


^travel 


3.3 X 10^^ 
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^travel 
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Table 2 — Continued 
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P 






L 
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f ■ ^ 

t-min 


minl'f ■ f^T.^*^ 
-'-'■***-'-\''mini''ep J 




Flvpnt Ratp 

J — i V cxxij J. IjOjU^ 






[G] 


[ms] 


[erg s 








[cm] 






[erg s~i] 


[km-2 yr-l] 


Model 


F5 


1012 


5 


1 


.5 


X 


1040 


102 


10"' 


8.9 X 10^''^ 


^travel 


^travel 


8.8 X 1028 


9.5 X 10"'^ 


Model 


F6 


1012 


5 


1 


.5 


X 


1040 


102 


10'' 


8.9 X 10^^^ 


^travel 


^travel 


8.8 X 102"^ 


3.6 X 10"'^ 


Model 


F7 


1012 


5 


1 


.5 


X 


1040 
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10^ 


8.9 X 10^^ 


^travel 


^travel 


8.8 X 1026 


6.2 X 10^7 
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F8 


1012 


5 


1, 


.5 


X 


1040 
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8.9 X IQi^ 


^travel 


^travel 


8.8 X 1025 


7.6 X 10-8 
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Gl 


1012 


5 


1, 


.5 


X 


1040 


IQi 


IQi 


8.9 X 10^3 


^travel 


^travel 


3.4 X 10^2 


2.1 X 10-1° 
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G2 
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5 


1, 


.5 


X 
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IQi 
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3.0 X 10^^ 


^travel 


^travel 


4.3 X 10^1 


1.0 X 10-'^ 


Model 


G3 


1012 


5 


1, 


.5 


X 


1040 


IQi 
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3.0 X lO^"^ 


^travel 


^travel 


2.1 X lO^o 


7.2 X 10-''' 
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G4 
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5 


1 


.5 


X 
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IQi 
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^travel 
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Model 


G5 
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5 


1 


.5 


X 
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IQi 
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^travel 


^travel 


2.1 X 1028 


2.3 X 10-'^ 
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G6 
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5 


1 


.5 


X 
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2.1 X 102^^ 


8.8 X 10^7 
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G7 
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5 


1, 


.5 


X 
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1.5 X 10-7 
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G8 
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5 


1, 


.5 


X 
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Shortest timescale among proton's synchrotron cooling timescale tsync , traveling timescale ttravel i ^^'^ collision timescale 
tco\- When tmin = ttravel or tool, it means that tmin is always shorter than tsync- On the other hand, when tmin = tcob it 
means that tool is shorter than any other timescales at some region of the nebula flow. 

^ Shorter timescale between tniin and energy transfer timescale from protons to electrons tep- When min(tinin,tep) = tmin, 
it means that tmin is always shorter than tep in the nebula flow. On the other hand, when min(tinin,tep) = tep, it means that 
tep is shorter than tmin at some region of the nebula flow. 

Note. — Luminosity of pions (ivr) and event rate at a km^ high-energy neutrino detector are not shown when the conditions 
(i)-(iii) presented in section 3 are not satisfied. 



